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SUMMARY 

Heterogeneity within pluripotent stem cell (PSC) populations is indicative of dynamic changes that occur when cells drift between 
different states. Although the role of metastability in PSCs is unclear, it appears to reflect heterogeneity in cell signaling. Using the Fucci 
cell-cycle indicator system, we show that elevated expression of developmental regulators in Gl is a major determinant of heterogeneity 
in human embryonic stem cells. Although signaling pathways remain active throughout the cell cycle, their contribution to heteroge- 
neous gene expression is restricted to Gl. Surprisingly, we identify dramatic changes in the levels of global 5-hydroxymethylcytosine, 
an unanticipated source of epigenetic heterogeneity that is tightly linked to cell-cycle progression and the expression of developmental 
regulators. When we evaluated gene expression in differentiating cells, we found that cell-cycle regulation of developmental regulators 
was maintained during lineage specification. Cell-cycle regulation of developmentally regulated transcription factors is therefore an 
inherent feature of the mechanisms underpinning differentiation. 



INTRODUCTION 

Pluripotent stem cells (PSCs) are heterogeneous under self- 
renewing conditions in culture (Enver et al., 2009; Graf and 
Stadtfeld, 2008; Martinez Arias and Brickman, 2011) and 
during embryonic development (Chazaud et al., 2006). 
This heterogeneity extends not only to the expression of 
pluripotency factors such as NANOG, REX1, and STELLA 
(Chambers et al., 2007; Hayashi et al., 2008; Singh et al., 
2007; Toyooka et al., 2008), but also to lineage-specific 
factors such as HEX, HES1, and GATA6 (Canham et al., 
2010; Kobayashi et al., 2009; Singh et al., 2007). Variations 
in gene expression are transient and reversible, indicating 
that PSCs alternate between different cell states. Although 
the function and molecular mechanisms underpinning 
this heterogeneity are unclear, it appears to be influenced 
by variations in the activity of signaling pathways at the 
single-cell level. WNT, BMP, NODAL, and FGF signaling 
through their downstream effectors has been implicated 
in contributing to PSC heterogeneity and serves to 
prime cells for differentiation when transiently activated 
(Galvin-Burgess et al., 2013; Price et al., 2013). As an 
example, heterogeneity can be significantly reduced 
when murine PSCs are cultured in the presence of small- 
molecule compounds that block ERK and GSK3 signaling 
(2i media) (Marks et al., 2012; Wray et al., 2011; Ying 
et al., 2008). In human embryonic stem cells (hESCs), sup- 
pression of WNT activity reduces signaling heterogeneities 



and the sporadic expression of developmental regulators 
such as BRACHYURY (Blauwkamp et al., 2012; Singh 
et al., 2012). Together, these observations indicate that 
signaling heterogeneities reflect alternate cell states that 
represent different differentiation potentialities. 

PSCs exhibit an unusual mode of cell-cycle regulation 
with a truncated Gl and a large percentage of S phase cells 
(Singh and Dalton, 2009). As PSCs differentiate, the cell 
cycle is remodeled, such that Gl is lengthened and the rela- 
tive amount of time associated with S phase cells is reduced. 
Recent reports (Calder et al., 2013; Coronado et al., 2013; 
Pauklin and Vallier, 2013) further documented this using 
the fluorescent ubiquitination-based cell-cycle indicator 
(Fucci) system (Sakaue-Sawano et al., 2008). Together, these 
studies point toward a direct relationship between the cell 
cycle and differentiation, consistent with earlier reports 
describing the ability of PSCs to initiate their differentia- 
tion program from Gl phase (Chetty et al., 2013; Jonk 
et al., 1992; Mummery et al., 1987; Sela et al., 2012; Singh 
and Dalton, 2009). This raises the possibility that heteroge- 
neous gene expression and cell signaling variations in PSCs 
may also be linked to cell-cycle progression. To address this 
question, we utilized the Fucci system in hESCs in combi- 
nation with fluorescence-activated cell sorting (FACS), 
and performed RNA sequencing (RNA-seq) analysis to 
establish that heterogeneous expression of developmental 
regulators is closely coupled to cell-cycle positioning. 
Our findings provide a rationale for gene-expression 
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Figure 1. Establishment and Character- 
ization of Fucci hESCs 

(A) Diagram of Fucci reporters. 

(B) Diagram of ceLL-cycLe phase-specific 
fluorescent protein expression patterns. 

(C) WA09 Fucci hESCs were imaged every 
15 min over a 24 hr period. Micron bar: 
50 u.m. Arrowheads: cell progression from 
GltoStoG2. 

(D and E) WA09 Fucci hESCs were treated 
with Violet Vybrant DyeCyde for 2 hr and 
then analyzed on the MoFLo. The approxi- 
mate percentage of cells in each cell-cycle 
phase is shown. Data are representative of 
biological replicates. 



DNA content 



heterogeneity in hESCs and a potential mechanism for 
"lineage priming" in Gl phase. Moreover, we show that 
transient activation of developmental genes in Gl, such 
as GATA6 and SOX17, is associated with the upregulation 
of 5-hydroxymethylation (5hmC), an epigenetic mark 
with proposed roles in gene activation (Szulwach et al., 
2011a; Yuet al., 2012b). 



RESULTS 

Establishment and Validation of Human Fucci ESCs 
To construct Fucci hESC lines, fluorescent reporters 
(Sakaue-Sawano et al., 2008) were introduced into expres- 
sion vectors under the control of the constitutive CAGi 
promoter, linked to either a neomycin (neo R ) or puro- 
mycin (puro R ) selectable marker through an internal 
ribosome entry site (Figure 1A). The principle behind 
Fucci indicators is that the "destruction box" sequences 
of the cell-cycle regulators CDT1 and GEMININ are fused 
to fluorescent proteins and, upon entry into a specific cell- 
cycle phase, are targeted for degradation. This allows for 
Fucci indicators to be present only during a specific 
window of time, with the advantage of reporting cell-cycle 
position in living cells. The expected pattern of reporter 
activity is summarized in Figure IB, showing an initial 
double-negative (DN) population indicative of early Gl 
cells, followed by accumulation of CDT1-K02 (Kusabira 
orange-2) reporter fluorescence (red) throughout the rest 
of Gl. As the cells enter S phase, K02 fluorescence is ex- 



tinguished and GEMININ-Azl (Azami green-1) fluores- 
cence increases until mitosis is completed. After drug 
selection and expansion of WA09 hESCs, we confirmed 
the authentic cell-cycle-regulated expression of Fucci 
indicators by using several approaches in live cells. Both 
fluorescent reporters were robustly expressed in a largely 
nonoverlapping pattern, consistent with their being 
associated with separate phases of the cell cycle (Figures 
IB and 1C). We next validated the phase-specific expres- 
sion of fluorescent reporters based on DNA content using 
Vybrant DyeCycle stain. Cells expressing the Kusabira 
orange-2 (K02) reporter were exclusively in Gl, whereas 
cells expressing the Azami green-1 (Azl) reporter were 
composed of S phase and G2/M cells. As the cells transi- 
tioned through S phase and G2/M, Azl reporter intensity 
increased (Figure ID). 

The total cell population stained with DyeCycle stain 
exhibited a typical hESC DNA profile comparable to that 
seen with other dyes, such as propidium iodide and 
Hoechst 33342 (Singh and Dalton, 2009), with a high per- 
centage (~43%) of S phase cells (Figure ID). Upon gating of 
the appropriate fractions, we confirmed the cell-cycle 
expression of fluorescent reporters. Specifically, the DN 
fraction represented early Gl cells, the K02 (red) fraction 
represented late Gl, the Azl-low (AzL, low green) repre- 
sented S phase, and the Azl -high (AzH, high green) repre- 
sented G2/M (Figure IE). Less than 5% of the cells were 
double positive for K02 and Azl, indicative of cells in early 
S phase. These data validate the Fucci reporters as faithful 
indicators of cell-cycle position, as judged by DNA content. 
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The Heterogeneity of hESCs Is Dependent on Cell- 
Cycle Position and Cell Signaling 

Although PSCs are known to have an unusual mode of cell- 
cycle regulation, deep mechanistic insight into this is lack- 
ing. Moreover, no global analysis of cell-cycle dependent 
genes has been performed in pluripotent cells. To examine 
this issue, we performed RNA-seq analysis on the four 
FACS-isolated Fucci cell-cycle fractions (DN, K02, AzL, 
and AzH). We identified ~500 transcripts that followed a 
reproducible pattern of cell-cycle regulation in three bio- 
logical replicate experiments (Figure 2 A; Table SI available 
online). Using unsupervised cluster analysis, we identified 
ten cohorts of transcripts, all of which displayed a similar 
pattern of periodicity, and then subjected them to Gene 
Ontology (GO) analysis (Figure 2B; Table S2). Surprisingly, 
the largest group of cell-cycle-dependent transcripts con- 
sisted of developmental regulators (Figure 2C). 

RNA-seq analysis indicated that pluripotency regulators 
such as OCT4, NANOG, and S0X2 showed no consistent 
pattern of periodicity in the cell cycle (Figures 2A and 
2D). Cell-cycle regulators associated with mitosis, such as 
CYCLIN El, CDK1, and AURORA A, showed a consistent 
pattern of cell-cycle regulation, but regulators of Gl and S 
phase progression, such as CYCLIN E and E2F1, lacked peri- 
odicity at the transcript and protein levels (Figures 2E, SI A, 
and SIB). To independently validate these findings, we 
sorted live hESCs based on DNA content using Vybrant 
DyeCycle (Figures S1C-S1E). Although it is not possible 
to separate early and late Gl cells using this approach, tran- 
script analysis of cell-cycle genes in Gl, S, and G2/M can 
still be reliably performed. Quantitative RT-PCR (qRT- 
PCR) analysis of Vybrant DyeCycle FACS-isolated hESCs 
generated results similar to those obtained with Fucci- 
sorted cells (Figure SI). This further validates the Fucci indi- 
cator hESC line as a reliable tool for studying the cell cycle. 

The most pronounced class of transcripts that showed 
cell-cycle-dependent periodicity was composed of develop- 
mentally regulated transcription factors. Transcripts identi- 
fied by GO analysis as being involved in development, 
cell-fate commitment, and transcriptional regulation 
were heavily enriched in cluster 6, coinciding with a cohort 
of genes with elevated expression in Gl phase (DN and 
K02) (Figure 2C). To confirm these findings, we assayed 
the transcript levels of several randomly selected develop- 
mental regulators from cluster 6 (Gl cohort) by qRT-PCR 
analysis. This analysis confirmed that transcripts impor- 
tant for embryonic patterning (HOXD1 and HOXD13), 
neural development (ZIC1, GBX2, EN1, PAX7, and 
SOX 11), and mesoderm/definitive endoderm (DE) specifi- 
cation (GATA4, TBX5, and PITX2), were all selectively upre- 
gulated in Gl phase (Figure 2F). These data show that 
hESCs in Gl have higher expression of transcription 
factors that are important for differentiation and develop- 



ment, indicating that hESCs exist in a "lineage-primed" 
state during the Gl phase of the cell cycle. 

The selective upregulation of transcripts for develop- 
mental genes in Gl was surprising, but it raised the 
possibility that such regulation could be related to previous 
observations of heterogeneity in PSCs (Enver et al., 2009; 
Graf and Stadtfeld, 2008; Martinez Arias and Brickman, 
2011). To establish the relationship between cell cycle 
and hESC heterogeneity, we characterized in greater detail 
the developmental genes that are strongly expressed as part 
of DE differentiation and weakly expressed during PSC self- 
renewal (Singh et al., 2012). When we assayed Fucci hESC 
fractions by qRT-PCR, we found that GATA6, SOX17, 
FOXA2, and GATA4 transcripts were all upregulated in Gl 
(Figure 3A). To establish that transcriptional control of 
developmental regulators is a determinant of their cell- 
cycle regulation, we pulse labeled Fucci cells with ethynyl 
uridine (EU) for 1 hr and evaluated the levels of newly syn- 
thesized transcripts by qRT-PCR (Figure S2). The levels of 
nascent transcripts for developmental regulators from clus- 
ter 6 were all clearly upregulated in Gl, consistent with 
global RNA-seq and qRT-PCR analyses (Figure 2F; Table 
SI). This indicates that transcriptional control of develop- 
mental genes is a major contributor to cell-cycle-dependent 
heterogeneity in hESCs. Immunostaining confirmed these 
results, demonstrating that developmentally regulated 
transcription factors are expressed during a narrow window 
of time in the hESC cell cycle (Figure 3B). Analysis at the 
single-cell level showed that developmental factors are 
coexpressed with pluripotency factors, indicating that 
spontaneously differentiating cells are not responsible for 
this pattern of heterogeneity (Figure S2). Expression of 
developmental regulators is extinguished as cells transition 
into S phase, thereby establishing a strong link between 
cell-cycle position and PSC heterogeneity. To further 
explore the cell-cycle regulation of developmental genes, 
we differentiated Fucci hESCs toward neuroectoderm pro- 
genitors using dual inhibition of ACTIVIN and BMP 
signaling (Chambers et al., 2009). After 5 days of differenti- 
ation, the expression of neural genes emerged in a cell- 
cycle-regulated pattern, peaking in early Gl (Figure S2). 
These data support the general idea that developmental 
regulators are cell-cycle regulated and contribute to hetero- 
geneity in PSCs, regardless of which germ layer these 
factors are associated with. 

Heterogeneity in PSCs has been largely attributed to 
stochastic changes in signaling networks, but up to now, 
no direct links to the cell cycle have been established. 
The data we obtained thus far introduced another factor 
that contributes to PSC heterogeneity: the cell cycle. 
Next, we sought to determine the relationship between 
cell signaling and cell-cycle position with respect to 
heterogeneous expression of developmental genes. We 
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Figure 2. Gene-Expression Profiling of Human Fucci ESCs 

(A) Following RNA-seq, cluster analysis was performed on cell-cycle-regulated transcripts in WA09 Fucci hESCs, represented by a heatmap. 

(B) GO analysis of ceU-cycle-regulated genes from human Fucci RNA-seq. 

(C) Percentage of genes by GO analysis in the transcriptional regulation or cell-cycle categories according to RNA-seq cluster analysis. 
Developmental transcription factors were heavily enriched in cluster 6, and mitotic genes were enriched in clusters 2 and 7. 

(D-F) WA09 Fucci hESCs were separated into cell-cycle phases by FACS, and RNA levels were then determined by qRT-PCR for (D) 
pluripotency genes, (E) cell-cycle mitotic genes, and (F) developmental transcription factors. 

Data are representative of biological triplicate experiments displayed as mean ± SEM. *p < 0.05, **p < 0.01. See also Figure SI. 



Stem Cell Reports | Vol. I | 532-544 | December 1 7, 20 1 3 ©20 1 3 The Authors 535 



Stem Cell Reports 

Cell-Cycle-Dependent Heterogeneity in hPSCs 



10 

c 

CD O 

> '« 

Ct- x 
LU 0 



GA7/46 



SOX? 7 



FOXA2 



GATA4 



ESC 
ESC+UO 
ESC+DKK1 




E-G1 L-G1 S G2 



E-G1 L-G1 S G2 



Az1 


K02 


NANOG 




0 

4 










Az1 


K02 

• 


FOXA2 

r 


» 


m 


> 






> 


Az1 


K02 

A 

% • 


soxw 

• * A 




V 

• • • 


V 




A 


A 



E-G1 L-G1 S G2 



,120 



E-G1 L-G1 S G2 



§ 100 
° 80 
° 60 
& 40 

CO 

£ 20 
0 



ESC ESC ESC 
+U0 +DKK1 



SOX17 
FOXA2 
NANOG 



WA09 hESCs 



E-G1 L-G1 


S 


G2 










NANOG 








pAKT S473 








pERK1,2 T202Y204 








pGSK3|3 S9 








GSK3p 








pS6 S235,236 

CDK2 



G1 



G1 



GATA6 fl 
SOX17/ I 
FOXA2 / I 



A_ 



o#ooo«o 



signal-responsive window (G1) 



Figure 3. WNT and ERK Signaling Promote Heterogeneity in Late Gl Cells 

(A) WA09 Fucci hESCs grown in self-renewal media were treated with U0126 (20 u,M) or DKK1 (150 ng/mL) for 3 days and then separated 
into cell-cycle phases by FACS. RNA expression analysis by qRT-PCR was performed for endoderm markers. All experiments were performed 
in technical triplicate and are representative of multiple experiments. 

(B and C) WA09 Fucci hESCs were treated with U0126 (20 u.M) or DKK1 (150 ng/ml) for 3 days and immunostaining was performed for 
indicated proteins and quantitated in (C) as the ratio of Alexa FLuor 647+ (for indicated proteins) to K02+ cells. Data represent analyses 

(legend continued on next page) 
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specifically addressed this by examining whether signaling 
pathways activate developmental genes in manner that 
depends on the cell-cycle-phase. Fucci hESCs were treated 
with MEK/ERK signaling inhibitor (U0126) or WNT antag- 
onist (DKK1) for 3 days in self-renewal media, and then the 
levels of developmental regulators were examined by qRT- 
PCR and immunostaining (Figures 3A-3C). Inhibiting 
either ERK or WNT signaling reduced transcript levels for 
the developmental regulators GATA6, SOX17, F0XA2, 
and GATA4 in K02+ (late Gl) cells. In addition, treatment 
with U0126 or DKK1 in self-renewal media also led to a 
significant reduction in background protein levels for 
FOXA2 and SOX17, but had no effect on NANOG. These 
data indicate that WNT and ERK signaling pathways pro- 
mote heterogeneity in hESCs by activating the transcrip- 
tion of developmental regulators in Gl phase. It should 
be noted that although 70%-80% of K02+ cells express 
SOX17 protein, we observed a small percentage of cells 
(<2% of total cells) that expressed SOX17 in all other frac- 
tions of the cell cycle. We believe these cells to be spontane- 
ously differentiated cells. When we evaluated the activity 
of several signaling pathways in Fucci fractions, we 
observed no general pattern of Gl -dependent cell signal- 
ing (Figure 3D). For example, we found no evidence of 
phospho-dependent regulation of AKT1, ERK1/2, or 
GSK33 at different cell-cycle positions. In contrast, phos- 
phorylation of ribosomal S6 protein was upregulated in S 
and G2/M phases, as described previously (Shah et al., 
2003), but not in Gl. 

The above findings rule out the possibility that Gl- 
associated heterogeneity is driven by cell-cycle-regulated 
activation of signaling pathways. Instead, Gl appears to 
represent a brief window of time in which cells are compe- 
tent to respond to external signals at the level of gene acti- 
vation (Figure 3E). How cells sense differentiation signals as 
a function of the cell cycle is presently unclear, but PSCs are 
known to be particularly susceptible to differentiation cues 
in Gl (Chetty et al., 2013; Mummery et al., 1987; Sela et al., 
2012). Using Fucci hESCs, we then went on to confirm 
these previous findings (Figure S3). Fucci hESCs were sorted 
into cell-cycle fractions and replated in serum-containing 
differentiation media. After 48 hr, qRT-PCR analysis was 
performed for differentiation markers from all three germ 
layers. Cells in early and late Gl displayed increased marker 



expression compared with those from S or G2/M phases. To 
further investigate whether Gl cells are more sensitive to 
differentiation signals over other phases of the cell cycle, 
we treated Fucci hESCs for 3 hr in endoderm differentiation 
conditions, followed by FACS and qRT-PCR (Figure S3). We 
found that the 3 hr treatment of Fucci hESCs in differenti- 
ation conditions was sufficient to increase the expression 
levels of endoderm genes in Gl cells, but only weakly 
increased gene-expression levels in S and G2. These data 
provide further evidence that cells preferentially initiate 
their differentiation from the Gl phase of the cell cycle. 

Cell-Cycle-Dependent Heterogeneity Is Retained 
during Early Differentiation 

Our data show that many developmental regulators are 
heterogeneously expressed in hESCs by a cell-cycle- 
dependent mechanism. It was not clear, however, whether 
these developmental regulators maintain a heterogeneous 
pattern of expression linked to the cell cycle during differ- 
entiation. To examine this question, we utilized the 
well-established DE differentiation model, in which phos- 
phatidylinositol 3 kinase (PI3-kinase) signaling is limited 
and ACTIVIN A/SMAD2,3 signaling is increased (D'Amour 
et al., 2005; McLean et al., 2007; Singh et al., 2012; Teo 
et al., 2011). Fucci hESCs were differentiated to DE for 2 
or 4 days, and then cell-cycle fractions were isolated by 
FACS and subjected to transcript analysis by qRT-PCR (Fig- 
ures 4A and 4B). By day 2 (d2), the DE transcripts GATA6, 
SOX17, and FOXA2 were significantly elevated in all cell- 
cycle fractions, but most noticeably in Gl (K02; Figure 4C). 
Corresponding increases in cell-cycle-regulated SOX1 7 and 
FOXA2 protein were also observed at this time (Figures 4D- 
4F). The initial stage of DE differentiation is therefore 
characterized by elevated expression of lineage markers 
(GATA6, SOX17, and FOXA2) while cell-cycle regulation is 
retained. BRACHYUR Y however, does not seem to be under 
this mode of control and shows no obvious cell-cycle regu- 
lation in hESCs (Figures 4C and S2B). By day 4 (d4), the 
endoderm transcript levels increased further, and although 
they were higher in Gl, they were now also significantly 
elevated in S and G2/M phases. The transition from d2 to 
d4 under these conditions was associated with significant 
remodeling of the cell cycle, where K02-positive cells 
increased from <10% to >40% (Figure 4B). To rule out the 



from five different fields of view, and the total number of cells counted per condition was >350. Arrowheads: K02+ cells that express 
endoderm markers in self-renewal media, but not in the presence of the MEK/ERK inhibitor (U0126). Data are representative of multiple 
experiments. Micron bar: 50 u.m. 

(D) Immunoblotting (20 ug/lane) for FACS-isolated Fucci fractions for the indicated signaling molecules. 

(E) Schematic diagram indicating that cells in Gl are most responsive to differentiation-inducing signals. Colors are representative of 
Fucci fractions (gray, early Gl; red, late Gl; yellow, Gl/S; green, S-G2/M). Circles depict cells progressing through the cell cycle, with solid 
circles representing cells that express developmental regulators. 

Data are representative of biological triplicate experiments displayed as mean ± SEM. *p < 0.05, **p < 0.01. See also Figure S2. 
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Figure 4. Developmental Regulators Retain Cell-Cycle Regulation during Differentiation 

(A) Schematic for examination of DE celL-cyde-regulated transcripts. Fucci hESCs were differentiated to DE for 2 or 4 days, cell-cycle 
fractions were isolated by FACS, and RNA expression analysis was performed. 

(B) FLow-cytometry plots of undifferentiated ESCs, 2-day DE, and 4-day DE. 

(C) RNA expression analysis of Fucci cell-cycle fractions after differentiation to DE. All experiments were performed in technical triplicate 
and are representative of multiple experiments. 

(D and E) Immunostaining upon differentiation of Fucci ESCs for 2 or 4 days to DE is quantitated in (E) as the percentage of positive cells in 
each cell-cycle fraction. Data are representative of multiple independent experiments. Micron bar: 50 urn. 

(F) Immunoblotting (20 pg/lane) for the indicated proteins from cell-cycle fractions in undifferentiated ESC or after 2 days of DE 
differentiation. 

Data are representative of biological triplicate experiments displayed as mean + SEM. *p < 0.05, **p < 0.01. See also Figures S3 and S4. 
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Figure 5. The Epigenetic Mark 5hmC Is 
soxu Cell-Cycle Regulated 

(A) Dot blot analysis reveals a Loss of 5hmC 
in S phase. Representative of assays per- 
formed in biological triplicate. 

(B) MS analysis for 5hmC performed in 
technical triplicate and representative of 
experiments in biological duplicate. 

SOX2 (C) 5hmC-capture qPCR for upstream re- 

gions of TSS. 

(D) Sequencing of 5hmC-enriched regions 
reveals enrichment for DhMRs at TSSs and 
TESs. The experiment was performed in 

G2 technical triplicate and is representative of 
duplicate experiments. 

(E) DhMRs are colocated with enhancer 
marks (H3K4mel and K3K9ac) and bivalent 

. G2 marks (H3K4me3 and H3K27me3). 

s (F) TAB-seq for promoter regions in Fucci 
° l-g 1 fractions, with diagram indicating primer 
D E " G1 location. 

Data are representative of biological tripli- 
cate experiments displayed as mean ± SEM. 
*p < 0.05, **p < 0.01. See also Figure S5. 

It 
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possibility that late Gl (K02+) cells are simply a differenti- 
ated cell population, we examined the self-renewal and 
differentiation potential of these cells in greater detail (Fig- 
ure S4). After sorting late Gl cells and replating them in 
self-renewal or differentiation media, we found that they 
reestablished an asynchronous cell-cycle profile within 
2 days. From these replated late Gl cells, we sorted the 
Fucci cell-cycle fractions and analyzed gene-expression 
markers by qRT-PCR. Importantly, the endoderm markers 
reestablished their cell-cycle-regulated gene-expression 
patterns during self-renewal and differentiation. Immuno- 
staining of replated late-Gl cells further confirmed these 
results and showed clear differences in endoderm marker 
expression between early- and late-Gl cells. Altogether, 
our data indicate that although an increased proportion 
of cells in Gl may contribute to elevated levels of develop- 
mental regulators, other factors involved in the ampli- 
fication mechanism appear to play a role. Nevertheless, 
cell-cycle-dependent heterogeneity is retained as hESCs 
differentiate toward DE, indicating that transition through 
Gl is linked to preferential gene activation. 

Epigenetic Changes at Developmentally Regulated 
Genes in Gl Phase 

Although signaling pathways are activated in a cell-cycle- 
independent manner, cells seem to be responsive only in 
Gl phase, indicating that uncharacterized mechanisms 



link signal transduction to target genes required for cell- 
fate commitment. One mechanism involving signal- 
regulated transcription factors and cyclin D-associated 
cyclin-dependent kinase activities was recently described 
(Pauklin and Vallier, 2013). We hypothesized, however, 
that factors such as chromatin structure and epigenetic 
status may also be a corequirement for transcriptional acti- 
vation by signaling effectors in Gl. Because cytosine 
methylation (5-methylcytosine [5mC]) and 5-hydroxyme- 
thylcytosine (5hmC) are epigenetic modifications that are 
known to be important in modulating transcriptional 
activity (Szulwach et al., 2011a; Yu et al., 2012b), and little 
is known about 5mC/5hmC dynamics in the cell cycle, we 
decided to examine these marks in Fucci cells. To that end, 
we subjected genomic DNA from the four FACS-isolated 
Fucci hESC cell-cycle fractions to dot blot analysis and 
mass spectrometry to evaluate global 5hmC levels. Surpris- 
ingly, both analytical approaches revealed a clear upregula- 
tion of 5hmC in late Gl followed by a sharp decline during 
S phase (Figures 5 A and 5B). These data indicate that cell- 
cycle-regulated 5hmC oscillations occur globally during 
stem cell self-renewal. Furthermore, since 5hmC peaks in 
late Gl and is almost completely lost in S phase, this 
indicates that widespread DNA demethylation occurs in 
S phase. 

To further explore this issue in relation to developmental 
genes, we performed 5hmC-capture followed by qPCR on 
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Figure 6. PSC Heterogeneity Is Contingent upon the Cell-Cycle 
Position 

As cells pass through Gl phase, they become sensitized to signaling 
pathways and exhibit increased Levels of developmental regulators 
("lineage primed"). During differentiation, the expression of 
deveLopmentally important transcription factors increases, but 
they maintain periodicity during the cell cycle. Expansion of the Gl 
phase is Likely to be a factor in this. Colors are representative of 
Fucci fractions (gray, early Gl; red, late Gl; yellow, Gl/S; green, 
S-G2/M). Circles depict cells progressing through the cell cycle, 
with solid gray circles representing cells that weakly express 
developmental regulators, and solid black circles indicating cells 
that strongly express developmental regulators. 

genomic regions flanking GATA6, SOX1 7, CDK1, and SOX2 
(Figures 5C and S5). 5hmC was significantly enriched in 
late Gl (K02+) cells compared with all other fractions for 
GATA6 and SOX17, whereas CDK1 showed a slight enrich- 
ment in G2/M (AzH + ) cells. SOX2, however, did not have 
any significant enrichment between fractions. Changes 
in 5hmC at the GATA6 and SOX17 loci therefore parallel 
those seen in transcriptional activity. To further explore 
cell-cycle-regulated 5hmC regions globally, we performed 
5hmC-capture-seq in biological duplicate. A total of 
16,774 differentially 5-hydroxymethylated regions 
(DhMRs) that changed during the cell cycle were identified 
(Table S3). Importantly these 5hmC marks were enriched 
for regions encompassing transcription start sites (TSSs) 
and exons (Figure 5D). We further compared how these 
DhMRs overlapped with epigenetic marks at enhancers 
(H3K4mel and H3K9Ac) and bivalent marks (H3K4me3 
and H3K27me3) previously identified in hESCs (Figure 5E). 
Compared with randomly selected regions, the DhMRs 
were enriched for these epigenetic marks, suggesting that 
they have broad roles in chromatin organization and links 
to cell-cycle progression. 

To examine 5hmC cell-cycle changes with basepair reso- 
lution, we performed Tet-assisted bisulfite sequencing 
(TAB-seq) (Yu et al., 2012b) at the GATA6 and SOX17 loci 
(Figure 5F; Table S3) over regions where peaks were identi- 



fied by 5hmC-capture qPCR (see Figure 5C). Importantly 
single-base resolution TAB-seq confirmed our findings 
from the 5hmC-capture qPCR, i.e., 5hmC was enriched 
in the K02 fractions and declined as cells entered S phase. 
Although the exact function of 5hmC is still unclear, these 
data demonstrate two important points: (1) epigenetic 
marks are cell-cycle regulated, and (2) cell-cycle regulation 
of 5hmC may have important roles in regulating gene 
expression or chromatin structure. 



DISCUSSION 

PSCs exist in a dynamic, metastable state characterized by 
the heterogeneous expression of transcriptional regulators. 
Although variations in cell signaling contribute to these 
variations, it was not previously clear how heterogeneity 
within a cell population was generated. In this report, 
we show that developmental regulators are transiently 
expressed by a cell-cycle-dependent mechanism in hESCs 
and that Gl represents a narrow window of time in which 
hESCs are responsive to differentiation signals (Figure 6). 
By contrast, we find no cell-cycle regulation for pluripo- 
tency factors in hESCs, indicating that the cell cycle may 
not regulate the heterogeneity of factors such as NANOG 
or REX1. We propose that under self-renewing conditions, 
differentiation signals are not of sufficient magnitude or 
duration to elevate developmental regulators to a point 
where their expression can support lineage specification. 
Instead, low levels of signaling can only transiently activate 
transcription in Gl. It follows, then, that when hESCs are 
exposed to conditions that efficiently promote differentia- 
tion, such as elevated ACTIVIN A and low PI3K signaling, 
the magnitude and duration of transcriptional activation 
will be sufficient to activate a sustained developmental 
program. 

Why does Gl represent a special window of time that 
allows for signaling pathways to be sensed and develop- 
mental genes to be activated in response? This question is 
difficult to answer, but one possibility is that changes in 
chromatin structure in Gl are compatible with the activa- 
tion of signal-regulated genes (see Dalton 2013; Pauklin 
and Vallier, 2013). Although this issue requires further 
attention, in this report we show that one potentially 
important level of epigenetic modification, 5hmC, is cell- 
cycle regulated and changes in a cell-cycle-regulated 
pattern that follows the expression of developmental regu- 
lators such as GATA6 and SOX17. The elevated levels of 
5hmC combined with increased levels of transcription for 
developmentally regulated genes suggest that chromatin 
may be in a permissive state in Gl, which allows it to be 
more responsive to specification cues. This is consistent 
with the idea that cells initiate programs of differentiation 
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from the Gl phase. Whether 5hmC is a requirement or a 
consequence of the transcriptional changes seen in Gl 
will be a subject of continued investigation. Since GATA6 
and SOX 17 transcripts displayed significant levels of 
periodicity during the cell cycle, we chose to focus on these 
genes for further epigenetic analysis. Additional experi- 
mentation will be required to determine whether cell-cycle 
regulation of 5hmC reflects a general occurrence for all 
cell-cycle-regulated developmental genes. 

Another surprising observation described in this report is 
that, rather than being constitutively expressed with 
respect to the cell cycle during differentiation, develop- 
mental regulators retained a strong expression signature 
linked to transition through Gl. This again supports the 
idea that Gl represents a temporal window in which cells 
preferentially respond to signals, allowing them to activate 
transcription of developmental genes more efficiently. 
During the first 2 days of DE differentiation, the expression 
of GATA6 and SOX17 increased >100-fold with respect to 
hESCs. This preceded the expansion of Gl phase, a charac- 
teristic of cell-cycle remodeling during PSC differentiation 
(Singh and Dalton, 2009; White and Dalton, 2005), and is 
likely to reflect strong, sustained signaling as a conse- 
quence of cells being switched to media optimized for 
differentiation. The simplest interpretation is that elevated 
thresholds of SMAD2,3 target gene binding more effi- 
ciently and activate developmental genes in Gl phase 
(Pauklin and Vallier, 2013; Singh et al., 2012). In asynchro- 
nous cultures, DE transcripts such as GATA6 and SOX17 
begin to increase by d2 and peak by d4. Although signifi- 
cant transcription of these genes already occurs at d2, 
this is not apparent in typical differentiation experiments 
because only 10% of cells are in Gl. 

Between d2 and d4, transcripts continue to increase, 
and over this period the proportion of cells in Gl 
increases dramatically, contributing to a spike in DE tran- 
scripts. It should not be overlooked that expression of 
developmental genes increases in all phases during differ- 
entiation, indicating that entry and exit from Gl repre- 
sent a favorable window of time for transcription rather 
than serving as an on/off switch. In contrast, expression 
of BRACHYURY is not linked to cell-cycle progression 
and has a different mode of regulation during differentia- 
tion and in hESCs. In addition, we find that some neural 
genes, such as PAX6 and SOX1, lack cell-cycle regulation 
in hESCs but impose a cell-cycle-regulated pattern in 
neural progenitors. Other neural genes, such as OTX2, 
are only weakly cell-cycle regulated. Although our data 
show that a large cohort of developmental regulators 
were cell-cycle regulated, a subset exhibited no obvious 
periodicity. 

In summary, our findings provide a clear explanation for 
heterogeneity in the expression of developmental regula- 



tors in hESCs. Moreover, we show that cells are more 
responsive to differentiation signals in Gl than in other 
phases of the cell cycle, and thus this can be considered 
to be a "lineage-primed" phase in which cells are more 
susceptible to differentiation signals. This model is consis- 
tent with other reports showing that pluripotent cells 
initiate differentiation from the Gl phase of the cell cycle 
(Chetty et al., 2013; Jonk et al., 1992; Mummery et al., 
1987; Pauklin and Vallier, 2013; Sela et al., 2012). 

EXPERIMENTAL PROCEDURES 

hESC Culture and Differentiation 

WA09 hESCs were grown as previously described in media contain- 
ing HEREGULIN (31 (10 ng/ml; Peprotech), ACTIVIN A (10 ng/ml; 
R&D Systems), and human LONGR3 insulin-like growth factor 1 
(IGF-1; 200 ng/ml; Sigma; Singh et al., 2012). Human recombinant 
DKK1 was supplied by R&D Systems and U0126 was supplied by 
LC Laboratories. WA09 hESCs were differentiated to endoderm 
by withdrawing HEREGULIN pi and IGF-1, and increasing the 
ACTIVIN A concentration to 100 ng/ml. Also, BIO (2 uM; 
Calbiochem) was used during the first 24 hr. WA09 hESCs were 
differentiated to all three germ layers by plating in differentiation 
media consisting of 20% fetal bovine serum (Sigma) in Dulbecco's 
modified Eagle's medium (DMEM)/F-12. Differentiation to 
neuroectoderm was performed as previously described (Chambers 
et al., 2009). WA09 Fucci hESCs were plated in 20% KnockOut 
Serum Replacement (Invitrogen) in DMEM/F-12 (Cellgro) with 
SB431542 (20 uM; Tocris) and murine Noggin (500 ng/ml; R&D 
Systems) for 5 days. 

RNA-Seq 

Human Fucci ESCs were sorted into early-Gl (DN), Iate-Gl 
(K02), S phase (AzL), and G2/M (AzH) fractions by FACS on a 
MoFlo (Beckman Coulter) in biological triplicate (three indepen- 
dent passages). RNA was extracted using TriZol (Invitrogen) and 
subjected to sequencing with an Illumina HiSeq instrument 
(Hudson Alpha). For each biological replicate, ~S0 M reads per 
sample were obtained by 2 x 50 nucleotide paired-end 
sequencing. All analyses were performed based on the human 
reference genome hgl9. Alignment of RNA-seq reads against 
known genomic annotations downloaded from the University 
of California Santa Cruz site (http://www.genome.ucsc.edu) was 
performed using Bowtie version 0.12.7 (Langmead et al., 2009) 
and TopHat version 1.3.3 (Trapnell et al., 2009). Transcript-level 
expression analysis was performed using Cufflinks version 1.2.1 
(Trapnell et al., 2010). Transcripts that showed consistent 
patterns of change in all three biological replicate experiments 
were selected for downstream analysis. Transcripts were filtered 
for Spearman correlation coefficients, with a cutoff of 0.5 
between any two of the three replicates. 

5hmC Analyses 

Dot blot analysis was performed as previously described (Yu et al., 
2012b). Briefly, genomic DNA was denatured with NaOH, 
neutralized with ammonium acetate pH 7.0. and spotted onto 
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a nitrocellulose membrane. Following blocking in 5% nonfat 
milk, filters were probed with 5hmC antibody (Active Motif). 
The filters were then probed with horseradish peroxidase-conju- 
gated secondary antibody and detected with enhanced chemitu- 
minescence reagents. Global 5hmC levels were assessed by liquid 
chromatography-electrospray ionization tandem mass spectrom- 
etry (LC-ESI-MS/MS) as previously described (Yu et al., 2012b). 
The amount of global 5hmC was expressed as the proportion 
of total cytosine. Briefly, 1 ug of genomic DNA was first dena- 
tured by heating at 100°C. Five units of Nuclease PI (Sigma) 
was added and the mixture was incubated at 37°C for 2 hr. A 
1/10 volume of 1 M ammonium bicarbonate and 0.002 U of 
venom phosphodiesterase 1 (Sigma) were added to the mixture 
and the incubation was continued for 2 hr at 37°C. Next, 
0.5 U of alkaline phosphatase (Roche) was added and the mixture 
was incubated for 1 hr at 37°C. Before injection into an Agilent 
Zorbax 4.6 mm x 50 mm, 3.5 um particle size -cl8 column, 
the reactions were diluted with water to dilute out the salts and 
enzymes. LC separation was performed at a flow rate of 220 
I/min. Quantification was done using a LC-ESI-MS/MS system 
in the multiple reaction monitoring (MRM) mode as previously 
described (Figueroa et al., 2010). 5hmC-capture-seq and TAB- 
seq were performed exactly as previously described (Szulwach 
et al., 2011b; Yu et al., 2012a). 

qRT-PCR, Immunoblotting, and Immunostaining 

All qRT-PCRs were performed using Taqman Assays on Demand 
(Applied Biosystems) on an iCycler (Bio-Rad). All assays were 
performed in triplicate, normalized to glyceraldehyde 3-phos- 
phate dehydrogenase, and analyzed using the AACT method. 
The data are representative of multiple experiments. Gene 
expression of newly synthesized transcripts was determined 
using the Click-iT Nascent RNA Capture Kit (Invitrogen) followed 
by qRT-PCR with Taqman assays. Fucci hESCs were treated for 
1 hr with EU and kept on ice before and after FACS. FACS 
isolation was performed at 4°C. Immunoblotting and immu- 
nostaining were performed as previously described (Singh 
et al., 2012) with antibodies against NANOG (GTX100863; 
GeneTex), OCT4 (sc-8628), pl07 (sc-318), pl30 (sc-317), CYCLIN 
A (sc-596), CYCLIN Bl (sc-752), CYCLIN E (sc-481), p27 KIP1 
(sc-527), Mem (sc-48407), CDK2 (sc-163; Santa Cruz Bio- 
technology), p-CDKl (9111S), pS6 (4858S), pAKTl (4060S), 
pERK (9101S), pGsk3|3 (9323S), GSK3fj (9315S; Cell Signaling 
Technology), SOX17 (AF1924; R&D systems), and FOXA2 
(07-633; Millipore). 

Time-Lapse Microscopy 

Fucci hESCs were imaged with an Olympus Viva View fluorescent 
microscope-incubator at 15 min intervals over a 24 hr period 
with 200 ms exposure time. 

Statistical Analyses 

All assays were performed in technical triplicate (unless otherwise 
indicated) and are representative of biological replicates. All data 
are displayed as mean ± SEM. Statistical analyses were performed 
using one-way ANOVA for all assays. 
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